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Abstract 

This paper presents elements about the radial orbit instability, which occurs in 
spherical self-gravitating systems with a strong anisotropy in the radial velocity 
direction. It contains an overview on the history of radial orbit instability. We 
also present the symplectic method we use to explore stability of equilibrium 
states, directly related to the dissipation induced instability mechanism well 
known in theoretical mechanics and plasma physics. 
Keywords : gravitation, instability, radial orbits, symplectic 

1 Introduction 

An interesting problem in Astrophysics is the study of A^-body self-gravitating 
systems with a lot of radial orbits, when most particles have very prolonged 
orbits with near-radial velocities, and come close to the system's center. The 
stability of such systems remains an open question to this day. although several 
elements suggest out that spherically symmetric systems of radial orbits could 
be unstable, and susceptible to lose their initial sphericity. This is called the 
radial orbit instability, henceforth referred to as Roi. This mechanism could 
actually be an explanation for the shape of some astrophysical objects, such as 
elliptical galaxies, that are hard to explain otherwise. 

For sufficiently large values of N, which is certainly the case for galaxies, A^- 
body self-gravitating systems can be described by a smooth distribution func- 
tion, and two-body interactions or "collisions" can be neglected in front of the 
whole gravitational potential created by this function^ For simplicity, we will 
assume that each body has the same mass m. 



* Electronic address : lionel.marechal@ensta.fr 
^Electronic address : jerome.perez@ensta.fr 

^To see that, we can consider the two characteristic time constants of its dynamics. Cross- 
ing time is the typical time a particle takes to go across the system during its movement, 
it describes its trajectory in the potential created by other particles as a whole. Relaxation 
time is the time it takes for a particle's trajectory to be significantly influenced by collisions. 
A calculation, based on the velocity change during an encounter and the typical number of 
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Combining the Vlasov equation (describing how a distribution function evol- 
ves in time in a given potential when collisions are neglected) and the Poisson 
equation (giving the potential created by this very function), we get the Vlasov- 
Poisson system : 

/ f + ^-Vq/-mVqV^.Vp/ = |£ + {/,ii;} = o 

\ V.(q) = -Gm /^dP' ^ > 

Here, dP is a shorthand for d'^qd'^p, and the notation a' does not denote a deriva- 
tive of a, but instead the value a(q', p', t) (for any variable a). E is the average 
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energy per particle E (q, p, t) := ^ + miplq) (or one-particle Hamiltonian) . 

Spherically symmetric, stationary solutions to this system can be either 
isotropic in velocity space, in which case the distribution function can be writ- 
ten as a function of the energy per particle alone, / = /o(-E); or, in the case 
where velocity distribution is not isotr opic, as a function of energy and an gular 
momentum fo{E, L"^) (see for example Perez. Alimi. Alv. and Scholll ll99 6lFl 



Considering the amount of work done about radial orbit instability (Roi) 
over the years, and the importance of this concept, we felt it was necessary to 
do a complete timeline of publications about it; this will be the second part of 
the present paper. In a third part, we will present the symplectic method that 
we use to study the stability of self-gravitating systems, and its prospects on 
the study of radial orbit instability. 



2 Historical overview of radial orbit instability 
2.1 The pioneers 

The fi rst important work published about Roi is an analytical result bv lAntonov 



(|1973[ ). It establishes a differential system for a given "displacement" of orbits 
and the corresponding Poisson equation, in the limit case of radial orbits. Insta- 
bility of the system is then proved by constructing a strict Lyapunov function 
for this system, al though the pr oof is unclear. 

The same vear. lHenon (I1973I) published one of the first numerical simulations 



of the problem, using TV = 1000 spherical shells. Isotropic, polytropic models 
f{E) oc E" were found to be stable (as it was obtained by Antonov in the 60s), 
whereas anisotropic systems, for generalized polytropes f{E) oc E^L^™, were 
found to become unstable when m — >■ — 1, which corresponds to a system with 
more and more radial orbits. The method does not allow to see the effect on the 
p osition spa ce (the famous "bar"). This article made no reference to the work 



of Antonov 



A contrary result was published by the French team of lGillon. Doremus. and Baumann 



( 19761 ). Using water bag methods (decomposing distribution functions as a sum 



encounters during one crossing, shows (see for instance iBinnev and Tremaind ll2008f) . chapter 
1.2.1.) that relaxation time is typically of order larger than crossing time. This means 

that for durations comparable to the crossing time, the effect of collisions can legitimately be 
neglected. 

^Conversely, a distribution function fo{E) describes indeed an isotropic, spherically sym- 
metric equilibrium, and fo{E,L^) a spherically symmetric, anisotropic equilibrium. 



2 



of functions constant over phase space domains) H they predicted that all self- 
gravitating systems f{E,L^) were stable against non spherical perturbations. 

2.2 Subsequent advances 

In the eighties, IPolvachenko and Shukhman (1981) proposed a matrix formu- 



lation of the stability problem, based on a Fourier series decomposition of 
perturbations; this allows them to prove that models of the form f{E — 

(later called Ossipkov-Merritt models) are unstable if is sufficiently small. 
The article unfortunately isn't very clear either, but the result is contrary to 
iGillon et al.l and confirms Antonov and Henon. This article also presents an 
often-quoted stability criterion, stating that radial orbit instability occurs when 
> 1.75 ± 0.25, Tr and Tj_ being respectively the system's total radial and 
orthogonal kinetic energies. 

The first full, rea listic numerical su rvey of the problem of gravitational col- 
lapse was realized bv lvan Albadal (1982). This survey considers sets of = 5000 



particles, for different initial conditions : homogeneous spheres and systems 
made of smaller homogeneous spheres (clumps) , for different speed distributions 
determined by the initial virial ratio. The results are clear : while collapsing 
homogeneous spheres are not affected by radial orbit instability, clumped sys- 
tems with violent collapse (small virial ratio) lead to a triaxial equilibrium. The 
resulting light and density profiles are compatible with the ones observed on 
galaxies when projected onto the sky : the de Vaucouleurs R^^^ where R is the 
projected radius. 

One of the first complete studies of instability, both analytical and numerical, 
was made by 'Barn es. Hut, and G oodmarJ (Il98d) . Numerical studies with A^- 



body methods confirm and complete Hcnon's results for generalized Plummer 
models, as well as the Russians' stability criterion. Their analytical explanation 
for Roi links it to Jeans instability, stellar pressure in the tangential direction 
would no longer be sufficient to offset the natural tendency of radial orbits to 
condense. 

The next paper on the subjec10 by Merritt and Aguilai ( 1985h . focuses on 



numerical results. It uses A^-body simulations with N = 5 ■ 10^, taking as 
initial conditions "galactic type" density profiles. Those models are Ossipkov- 
Merritt models (which are isotropic near the center and anisotropic towards the 
borders) generated from the well-known Jaffe model p{r) oc (?'/ro)^^(l-l-r/ro)^^ 
(compatible with the projected R^^^ profile). The results are that transition 
between stability and instability is fairly sharp, and happens for « 2.5, a 
bit more than predicted by the Russian criterion. However, comparison with 
a distribution function decreasing in E and in seems to point out that the 
value of isn't a reliable stability criterion. The article also puts forward 
(apparently for the first time) the idea that Roi can be useful in explaining 
galactic formation. 

Linear stability of isotropic spheres agains t all perturbations, an d of anisotropic 
spheres against radial ones was achieved by Kandrup and Svgnetl (fT98 5) using 



^Five years earlier, Doremus, Feix and Baumann had obtained the stabiUty of isotropic 
systems by the same met hod. 

''it actually references I Barnes et al.l lll986t) . in spite of an apparently earlier publication 
date. 
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energy methods based on an old principle by iKulsrud and Mark ()l970l) : this 
powerful and simple method requires the construction of a Hermitian operator. 
This technique seems irrelevant for the study of radial orbit instability, which 
is naturally associated to an anisotropic system receiving non-radial linear per- 
turbations, for in their case the operator is no longer Hermitian. However, a 
global or energetic ap proach of thi s problem was proposed by a short original 
but nebulous paper bv lMav and BinnevI |l986h . When it is applied to study the 
stability against non-radi a l p erturbations of a family of anisotropic isochrone 



spheres, iMav and Binnev (fl986') claim that it is compatible with the Russian 
criterion by Poly achenko and Shukhman. Two year later, a technical paper by 
Goodmanl (1988) clarifies the method proposed by May and Binney, but no 



application to the radi al orbit instability is proposed. 

The English team of lPalmer and Papaloizo ul ()l987l) made a purely analytica 



study (t he first one s ince the alleged Russian result, disregarding the waterbag 
result of lGillon et al. ) based on a spectral analysis of perturbations, decompos- 
ing them once again on a family of orthogonal functions. This study seems 
to indicate an instability, although is seems very hard to verify. Two other 
important aspects of this paper are a "demonstration" that the Russian crite- 
rion (already shaken by Merritt and Aguilar) is invalid, a nd a explana t ion of a 
mechanism for instability growth, inspired by a work by iLvnden-Belll (|l979l) : 
an bar-like perturbation of the potential in a spiral galaxy could influence a 
star's orbit and lengthen it, which tends to align orbits along the perturbation. 
This effect could play a part in bar formation in spir al galaxie s. 

A clear synthesis of all those results was made bv iMerritt (1987), including 
Lynden-Bell's mechanism, as well as a criticism of Jeans instability mechanism 
for needing an homogeneous system which is not the case here. 

iKatj (|l99ll ) brought in a new kind of simulations in this context, which 
showed that "authentic" cosmological simulations with a Hubble flow and merg- 
ing phases te nd to erase traces of possible primordial Roi. The same year, 
an article by ISahal (|l99ll) extended the reach of spectral methods for normal 
modes to infinite-extension mo dels, which was n ot the case for previous studies. 
Still the same year, a study bv IWeinberd (1991) used the matrix methods initi- 
ated by the Russian school of Polyachenko, and found some results again. This 
analysis was later the subject of a complete article by the Argentine team of 
Cincotta. Nunez, and Muzziol(|l996l) . who studied transformation of loop orbits 
into bow orbits , in the s tyle o f Lynden-Bell's mechanism and in accordance to 
the intuition of iMerrittI (|l987l ). 



2.3 A renewed interest 



Taking advantage of some of their analytical results, iPerez et al.l ([19960 P^'O" 
posed and tested a new stability criterion for self-gravitating systems based on 
the nature of the perturbations it is submitted to. This criterion is validated 
on Ossipkov-Merritt models applied to polytropes, the number of particles in- 
volved reaching for the first time reasonable values of iV « lO'* for the whole 
of simulations. In their analytical results, they explain how waterbag methods 
are lacki ng in the field of ra dial orbits, which may explain the now-discredited 



result of lGillon et a!] (Il976l). 



The German team of Theis and Spurzem ( 19991 ) undertook an extensive nu- 
merical study of Roi, using dedicated "Grape" machines for collapses of Plum- 
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mer spheres with varying initial temperature. Growth rate of Roi is largely 
affected by potential softening, and very little by variations of the number of 
particles. Those simulations highlighted a very long-term evolution (its time 
scale is relaxation time) of triaxial systems produced by Roi towards a more 
or less spherical system, an evolution that according to the authors is due to 
collisions. 

A systematic study of gravitati onal collapse wi t h test s for several numerical 
parameters (TV, softening, . . . ) bv lRov and Pere3 ( 2004f ) allowed, among other 
results, to stress that Roi is dependent on the presence of robust inhomogeneities 
in the precollapse system. Only several-scale collapses can lead to Roi : collapses 
of homogeneous spheres fail to do so. Those results were completed and refined 
bv lBoilv and Athanassoulal ( 20061 ). who showed a small effect of particle number 
on the final state. 

Although t he role of Ro i in structure formation was hurt by the aforemen- 
tioned work of Katz lll991 ), com plementary analyses by the German team of 



luss. Jain, and Steinmeta l 19991 ) and by the Canadian team of lMacMillan. Widrow. and Henriksen 



( 20061 ) observed the result of medium-scale structure formation by collapse ex- 
periments, with the possibility of numerically suppressing Roi. Acceptable den- 
sity profiles are only found when Roi actually takes place during primordial 
phases; otherwise structure profiles are incompatible with simulations and ob- 
servations! This gives a new argument for Roi as a fundamental process of 
structure formation. 

A new activity in this dom ain is ri sing since 2005, notably from E. Barnes's 
tea m. Their articles, notably Barnes. Williams, Babul, and Dalcanton (200 
and lBellovarv. Dalcanton. Babul. Quinn. Maas. Austin. Williams, and Barnei 
show that Roi not only creates a triaxial system in position space (which 
was known for a long time), but also creates a spatial segregation in veloc- 
ity space (with an isotropic center and a radial halo). This segregation could 
be the root of th e universal profile observed i n large structures. Follow-up 
papers, including iBarnes. Lanzel. and Williamd ([2009), indicated that Roi is 
not fou nd in constant-density c ollapse, a statement backed up by the Italian 
team of Tre nti and BertinI ( 20061) . This agrees with the previous explanation by 
Roy and Pe rez (200^ that Roi does not happen without a primordial equilib- 




rium state caused by a several-scale inhomogeneous collapse- 



La stly, Roi has been observed in a triaxial state by Antonini, Capuzzo-Dolcetta. and Merrittj 
( 20091 ) : it apparently happens when this state is populated with too many "box 
orbits with predominantly radial motions". The system would then become 
more prolate, and still triaxial. 

As this overview shows, Roi is a phenomenon that has been known for almost 
40 years, but still sparks some controversy and contradictory results. It seems 
to be fundamental in structure formation, yet there are still a lot of unanswered 
questions on both the physical and the analytical sense. Under which conditions 
does it happen, if it happens? Can we explain its mechanism? Several tools 
can be used to study radial orbit systems, and among them, potential energy 
methods could prove quite useful. We chose to focus on the symplectic method, 
which will be covered in the next section. 
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3 The symplectic method 



3.1 Presentation 



The symplectic method is a way of investigating the sta bihty of a steady stat e 
against possible perturbations. It was first developed bv [Bartholomew (1971), 
though its use in the study of gravita t ional pl asmas is a relatively recent de- 



Morrison and Kotschenreuther 



veloranent. See for instanc e iMorrisonI (jigs'? 
(I199(TI) and lKandrupl (Il99lh . 

This method makes use of the Hamiltonian structure of the system under 
scrutiny. The Vlasov-Poisson system indeed derives from the following Hamil- 
tonian : 



H [./] = j dr|^/ (r, t) - icm^ 1 dr 1 dr 



|q-q'l 



which is just the system's total energy. Note that its functional derivative ^ 

is the one-particle Hamiltonian E = ^ + mip. 

This structure allows us to easily compute the time variation of any func- 
tional of the distribution function. If K[f] is a derivable functional of /, then 
by definition of ^ : 



dK[f] 
dt 



SK df 



= TTi^dr= ^{i?,/}dr 



5f dt 



SK 



Sf 



(2) 



Here, we use the noncanonical Poisson bracket, which was introduced bv lMorrison and Greene! 
( 19801 ) for studies in hydrodynamics and magnetohydrodynamics (and thus 
sometimes referred to as the Morrison bracket). For two functional A and 
B of /, it is defined by 



Then, from (0) it follows that 



dt 



//{f.f}-^>'-H/. 



(3) 



As shown by Kandrup ( 199lh . all physical perturbations Z'-^' that /o can 
receive may be written in the form : 

/W iT,t) = -{gJo} 

g is called the generator of the perturbation. We denote as G the following 
operator : 

G[f] ■■= I fgdT 

Writing perturbations in this form allows us to compute the energy variation. 
A first-order calculation shows that 



i7W[/o] = [G,i/] (/o) 



5{/o,^;}dr = o 
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Since /o is a steady state, {/o, E} = so the energy variation is zero at first 
order. This corresponds to the classic definition of an equilibrium. 
The second-order variation is 

i/(2)[/o] = [G, [G,H]]{fo) 

The complete calculation of ij'^' is then possible. 



[G,i?](/)= I f{g,E}dT 



>From the expression of E, we get directly 
It leads to 



S[G,H] 
Sf 



= {9,E} 



__Grn± 

|q-q' 
Gm2 



dr' 



|q-q 
!^{.9'./'}dr' 



(4) 



>From this, computing i/*^^-' is easy 



H^'Hfo] = -J 



S[G,H] 



{5,/o}dr 



{9, E}{g, /o}dr - Gm^ J j ^''{^^''^['^ 



dVdV 



(5) 



This method is a very powerful one, since it allows us to obtain 77^^) effi- 
ciently in a very general case. Most methods to compute energy variation due 
to a perturbation require the knowledge of the system's actual geometry, that 
is to say /o(r, v) and i!(r)^ which are difficult to obtain knowing only fo{E, L^) 
(except in some particular cases). With this method, we do not need to know 
the system's geometry to compute the energy variation. 

3.2 Stability criterion 

In the classic case of one particle influenced only by conservative forces, the 
stability of an equilibrium is directly linked to the sign of the second-order 
energy variation : if it is positive semi-definite, then the equilibrium is stable; 
otherwise it is unstable. In the previous part, we derived the second-order 
energy variation around an equilibrium state : its sign should p rovide a cr i terion 
to know whether it is stable o r unst able (as seen, again, in iMorrisonL Il987 : 
Morrison and Kotschenreuther . 1990l) . Unfortunately, in more general cases, 
this is not so simpleH 

In the case where JJ *-^-* > for all generators g, then there is a definitive result 
bv [Bartholomew (1971), which proves that the system is stable. The symplectic 



^An example is the case of a charged particle in a negative harmonic potential, with a 
strong enough magnetic fi eld. This example is covered in more details in another paper 



jMarechal and Perezj . |2010D . 
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method at least allows to prove the stability of an equilibrium. However, in the 
case where there are generators g such that H'^'^^ < 0, there is no definite proof 

of instability, at least not without more hypotheses; 

An important mathematical result was given bv lBloch. Krishnaprasad. Marsden. and Ratid 
m a quite general case. Consider a Hamiltonian system with finite di- 



mension, that is initially at equilibrium . Then we suppose there exists a negative 
energy mode. In this case. lBloch et al.l proved that with the addition of dissipa- 
tion, the equilibrium becomes spectrally unstable, from which follow linear and 
nonlinear instability. This kind of instability can be called a dissipation-induced 
instability. 

As the Vlasov-Poisson is infinite-dimensi onal, this result doesn ' t dire ctly 
apply to our problem. More recent works by iRouslan and Marsden (I2OO9I) on 



the infinite dimensional case seem to indicate that it works in the same way; 
there is no general proof yet, but a result seems likely in the near future. 

The previous method allows us to retrieve previous results much more easily. 
In the isotropic case fo{E), this method proves stability against all perturba- 
tions provided that c^b/o < 0. In the anisotropic case fo{E,L'^), it proves (with 
the same condition OeJo < ) stability against all so- called preserving pert ur- 



bations, i.e. perturbations that verify {g,L^} = 0, see IPerez and Alvl ()l996l) 



In this context, rigorous and complementary results were proposed by a 



Germ an team with very different methods, and are compiled in Rein and Guol 



(|2003D . 



3.3 Perspectives for radial orbits 

A system with nearly radial orbits can be easily described with the previous 
formalism, as fo{E,L'^) = ip{E)5{L'^) with 5 a function that selects values near 
(for instance a Dirac distribution). In this case, it is possible to show that for 
a sufBciently selective function, there ar e negative energy modes. T he details 
have been published in another article (jMarechal and PerezL l201(l published 
online) . 

The idea is that H^'^'^ has two main components : 

= -j {.9, E}{g^ /}dr -Gm^ j j i^l/H^drdF' 

^ V V ' 

i.A) (B) 

The term [A) mostly corresponds to kinetic energy variation, while [B) corre- 
sponds to potential energy variation. For a distribution function that is suffi- 
ciently radial, it is possible to show that there is a class of generators g breaking 
spherical symmetry, such that {A) is negligible in front of {B). The latter is 
negative, as it can readily be written as the integral of a function against its 
own Laplacian. This proves the existence of negative energy modes. 



4 Conclusion 

We have shown a criterion to discuss the stability of Hamiltonian systems, and 
used it in the case of self-gravitating systems such as star clusters or galaxies, 
thanks to the the fact that Vlasov-Poisson is a Hamiltonian system. This cri- 
terion works on initial equilibrium states. The addition of dissipation seems 
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to imply an instability if a steady state has negative energy modes. Appar- 
ently, systems populated with radial orbits have negative energy modes, which 
would mean they are unstable and susceptible to lose their spherical symmetry, 
triggering radial orbit instability. 

It may be possible to give a physical interpretation of this result. Radial 
orbits have no tangential velocity, so they do not precess around the system 
center; instead they are confined on a line. Bringing the orbits closer is be 
possible, as they do not precess, and it leads to a lower energy state as the 
stars would be on average closer to each other than before. If a direction had a 
higher than average density, other orbits would tend to align in this direction, 
in a lower energy state; those orbits would bring on others, and so on, leading to 
an instability. This process would not be possible if the stars cannot dissipate 
energy; hence the necessity of dissipation to insure that radial orbit instability 
takes place. 



References 

F. Antonini, R. Capuzzo-Dolcetta, and D. Merritt. A counterpart to the radial- 
orbit instability in triaxial stellar systems. Mon. Not. R. Astr. Soc, pages 
1180-+, August 2009. doi: 10.1111/j. 1365-2966.2009. 15342.x. 

V. A. Antonov. On the instability of stationary spherical models with purely 
radial motio ns. In Dynamics of Galaxies and Star Clusters^ pages 139-143, 
translated in de Zeeuw (Il987l) . 1973. 



E. I. Barnes, L. L. R. Williams, A. Babul, and J. J. Dalcanton. Scale Lengths in 
Dark Matter Halos. ApJ, 634:775-783, November 2005. doi: 10.1086/497066. 

E. I. Barnes, P. A. Lanzel, and L. L. R. Williams. The Radial Orbit Instability 
in CoUisionless N-Body Simulations. ApJ, 704:372-384, October 2009. doi: 
10.1088/0004-637X/704/1/372. 

J. Barnes, P. Hut, and J. Goodman. Dynamical instabilities in spherical stellar 
systems. ApJ, 300:112-131, January 1986. doi: 10.1086/163786. 

P. Bartholomew. On the theory of stability of galaxies. Mon. Not. R. Astr. 
Soc, 151:333-+, 1971. 

J. M. Bellovary, J. J. Dalcanton, A. Babul, T. R. Quinn, R. W. Maas, C. G. 
Austin, L. L. R. Williams, and E. I. Barnes. The Role of the Radial Orbit 
Instability in Dark Matter Halo Formation and Structure. ApJ, 685:739-751, 
October 2008. doi: 10.1086/591120. 

J. Binney and S. Tremaine. Galactic Dynamics. Princeton University Press, 
2008. 

A.M. Bloch, P.S. Krishnaprasad, J.E. Marsden, and T.S. Ratiu. Dissipation 
induced instabilities. Ann. Inst. Henri Poincare, 1994. 

C. M. Boily and E. Athanassoula. On the equilibrium morphology of systems 
drawn from spherical collapse experiments. Mon. Not. R. Astr. Soc, 369: 
608-624, June 2006. doi: 10.1111/j. 1365-2966. 2006.10365.x. 



9 



p. M. Cincotta, J. A. Nunez, and J. C. Muzzio. On the Radial Orbit Instability. 

ApJ, 456:274-+, January 1996. doi: 10.1086/176647. 

P. T. de Zeeuw, editor. Structure and dynamics of elliptical galaxies; Proceedings 
of the lA U Symposium, Institute for Advanced Study, Princeton, NJ, May 27- 
31, 1986, volume 127 of lAU Symposium, 1987. 

D. Gillon, J. P. Doremus, and G. Baumann. Stability of self-gravitating systems 
with phase space density - A function of energy and angular momentum for 
aspherical modes. Astron. & Astroph., 48:467 474. May 1976. 

J. Goodman. An instability test for nonrotating galaxies. ApJ, 329:612-617, 
June 1988. doi: 10.1086/166407. 

A. Huss, B. Jain, and M. Steinmetz. How Universal Are the Density Profiles of 
Dark Halos? ApJ, 517:64-69, May 1999. doi: 10.1086/307161. 

Michel Henon. Numerical Experiments on the Stability of Spherical Stellar 
Systems. Astron. & Astroph., 24:229-+, April 1973. 

H. E. Kandrup. A stabiHty criterion for any coUisionless stellar equilibrium and 
some concrete applications thereof. ApJ, 370:312-317, March 1991. 

H. E. Kandrup and J. F. Sygnet. A simple proof of dynamical stability for a class 
of spherical clusters. ApJ, 298:27-33, November 1985. doi: 10.1086/163586. 

N. Katz. Dissipationless collapse in an expanding universe. ApJ, 368:325-336, 
February 1991. doi: 10.1086/169696. 

R. M. Kulsrud and J. W.-K. Mark. Collective Instabilities and Waves for Inho- 
mogeneous Stellar Systems. I. The Necessary and Sufficient Energy Principle. 
ApJ, 160:471-+, May 1970. doi: 10.1086/150448. 

D. Lynden-Bell. On a mechanism that structures galaxies. Mon. Not. R. Astr. 
Soc, 187:101-107, April 1979. 

J. D. MacMillan, L. M. Widrow, and R. N. Henriksen. On Universal Halos and 
the Radial Orbit Instability. ApJ, 653:43-52, December 2006. doi: 10.1086/ 
508602. 

L. Marechal and J. Perez. Radial orbit instability as a dissipation-induced 
phenomenon. Mon. Not. R. Astr. Soc, 405:2785-2790,2010. doi: 10.1111/j. 
1365- 2966.2010.16663.x. published online. 

A. May and J. Binney. Testing the stability of stellar systems. Mon. Not. R. 
Astr. Soc, 221:13P-17P, July 1986. 

D. Merritt. Stability of elliptical galaxies - Numerical experiments. In P. T. de 
Zeeuw, editor, Structure and Dynamics of Elliptical Galaxies, volume 127 of 
I AU Symposium, pages 315-327, 1987. 

D. Merritt and L. A. Aguilar. A numerical study of the stability of spherical 
galaxies. Mon. Not. R. Astr. Soc, 217:787-804, December 1985. 



10 



p. J. Morrison. Variational Principle and Stability of Nonmonotonic Vlasov- 
Poisson Equilibria. Zeitschrift fur Naturforschung, 42a: 11 15-1 123, 1987. doi: 
10.1086/166407. 

P. J. Morrison and M. Kotschenreuther. The Free Energy Principle, Negative 
Energy Modes and Stability. In V. G. Baryakhtar, V. M. Cliernousenko, N. S. 
Erokhin, A. B. Sitenko, and V. E. Zakharov, editor, IV International Work- 
shop on Nonlinear and Turbulent Processes in Physics, volume 2 of Nonlinear 
World, pages 910-932, 1990. 

P. J. Morrison and J. M. Greene. Noncanonical Hamiltonian density formulation 

of hydrodynamics and ideal magnetohydrodynamics. Physical Review Letters, 
45:790-794, September 1980. doi: 10.1 103/PhysRevLett.45. 790. 

P. L. Palmer and J. Papaloizou. Instability in spherical stellar systems. Mon. 
Not. R. Astr. Soc, 224:1043-1053, February 1987. 

J. Perez and J.-J. Aly. Stability of spherical stellar systems - I. Analytical 
results. Mon. Not. R. Astr. Soc, 280:689-699, June 1996. 

J. Perez, J.-M. Alimi, J.-J. Aly, and H. Scholl. Stability of spherical stellar 
systems - II. Numerical results. Mon. Not. R. Astr. Soc, 280:700-710, June 
1996. 

V. L. Polyachenko and I. G. Shukhman. General Models of Collisionless Spher- 
ically Symmetric Stellar Systems - a Stability Analysis. Soviet Astronomy, 
25:533-+, October 1981. 

G. Rein and Y. Guo. Stable models of elliptical galaxies. Mon. Not. R. Astr. 
Soc, 344:1296-1306, October 2003. doi: 10.1046/j. 1365-8711. 2003.06920.x. 

K. Rouslan and J. Marsden. Dissipation-induced instability phenomena in 
infinite-dimensional systems. Arch. Rational Mech. Anal., 194(2):611-668, 
November 2009. 

F. Roy and J. Perez. Dissipationless collapse of a set of N massive particles. 
Mon. Not. R. Astr. Soc, 348:62-72, February 2004. doi: 10.1111/j. 1365- 2966. 

2004.07294.x. 

P. Saha. Unstable modes of a spherical stellar system. Mon. Not. R. Astr. Soc, 
248:494-502, February 1991. 

C. Theis and R. Spurzem. On the evolution of shape in N-body simulations. 
Astron. & Astroph., 341:361-370, January 1999. 

M. Trenti and G. Bertin. Partial Suppression of the Radial Orbit Instability in 
Stellar Systems. ApJ, 637:717-726, February 2006. doi: 10.1086/498637. 

T. S. van Albada. Dissipationless galaxy formation and the R to the 1/4-power 
law. Mon. Not. R. Astr. Soc, 201:939-955, December 1982. 

M. D. Weinberg. A search for instability in two families of spherical stellar 
models. ApJ, 368:66-78, February 1991. doi: 10.1086/169671. 



11 



